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Abstract 

The set theory relations G, \, A, fl, and U have corollaries in sub- 
space relations. Geometric Algebra is introduced as the ideal framework 
to explore these subspace operations. The relations G, \, and A are easily 
subsumed by Geometric Algebra for Euclidean metrics. A short compu- 
tation shows that the meet (n) and join (u) are resolved in a projection 
operator representation with the aid of one additional product beyond the 
standard Geometric Algebra products. The result is that the join can be 
computed even when the subspaces have a common factor, and the meet 
can be computed without knowing the join. All of the operations can be 
defined and computed in any signature (including degenerate signatures) 
by transforming the problem to an analogous problem in a different al- 
gebra through a transformation induced by a linear invertible function 
(a LIFT to a different algebra). The new results, as well as the tech- 
niques by which we reach them, add to the tools available for subspace 
computations. 

1 Introduction 

Operations on subspaces are useful in applications everywhere. Geometric Al- 
gebra, the intriguing algebra promoted by David Hestenes to unify and simplify 
many areas of mathematics ^ ^, is introduced as the ideal framework to ex- 
plore subspace operations. A large repertoire of operations to compute subspace 
operations are made available by Geometric Algebra, but a few holes remain, 
notably with respect to the meet and join of subspaces. This paper should re- 
solve the outstanding issues. A section on preliminaries makes this treatment 
reasonably self-contained. Four subspace operations are introduced on equal 
footing, motivated by the set theory relations (\, A, H, and U), common usage, 
and applied needs. An overview of the four subspace operations for all signa- 
tures indicates that they are all fundamental and interconnected. Experts will 
note that Geometric Algebra deals with oriented subspaces and that in general 
there is no oriented solution to the meet and join problem|^. This problem 
is avoided by representing unoriented subspaces by projection operators. This 
allows us to extend the meet and join defined in the previous literature 
to give meaningful (nonzero) results for any subspaces. The price we pay for 
this extension is that the meet and join presented in this paper are not linear. 
This is not a disadvantage, because our meet and join agree with the previous 
literature except when the previous literature results are zero. Many operations 
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in this paper have an arbitrary scale and orientation. This is addressed in the 
penultimate section, where a geometrical significance can be given to the linear 
result of zero from the previous literature. 

2 Preliminaries 

This section includes the definitions and motivation for four subspace operations. 
Following the new definitions is a review of Geometric Algebra. 

2.1 The Subspace Operations 

The set theory operations of \, A, n, and U can be apphed to subspaces of 7?.". 
However, in general the outcome will not be a subspace, since the set theory 
operations do not respect the linear structure of the subspaces. Four subspaces 
operations are defined that are motivated by the set theory operations, but that 
respect the linear structure and thus always produce subspaces. Let A and B 
be two subspaces of TZ". 

1. The Meet Operation 

The meet of A and B is the set {x g TZ^ : x ^ A and x e S}. In words it 
is the largest common subspace. It shall be denoted by ^ fl S. 

2. The Join Operation 

The join of A and B is the set {x G 7^" : 3a. E A and 3h E B such that 
X — a + b}. In words it is the span of the two subspaces (i.e. the smallest 
common superspace). It shall be denoted hy AU B. 

3. The Difference Operation 

The difference of A and B is the set {a G ^ : Vb G a - b = 0}, provided 
S is a subspace of A. In words it is the orthogonal complement of B in A. 
It shall be denoted by A\B. 

4. The Symmetric Difference Operation 

The symmetric difference of A and B is the set {x G 7?." : 3a G ^ and 
3b G S such that x = a + b and Vc G .4 fl S x • c = 0}. In words it is 
the orthogonal complement of the meet in the join. It shall be denoted by 
A/\B. Clearly, Al^B = {AU B)\{A n B). 

This paper always uses the symbols \, A, n, and U to denote subspace 
operations. 

2.2 A Review of Geometric Algebra 

Geometric Algebra is based on Clifford Algebra over a real vector space. A Clif- 
ford Algebra is an algebra generated by the scalars and the elements of a vector 
space with a metric. The algebra has a linear, associative, distributive product, 
and the square of a vector is the squared length determined by the metric (for 
more details on Chfford Algebras see|^). Let W^i^^ be a (p-l-g-l-r)-dimensional 
real vector space with a set of linearly independent, mutually orthogonal vec- 
tors, {pi, ...,pp,qi, ...,q,,ri, ...,rr}, such that = 1, = -1, and = 0. 
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The Clifford Algebra and the Geometric Algebra generated by TZpi-^ shall be 
denoted TZp^q^r- 

The essential difference between a Geometric Algebra and a Clifford Algebra 
is that the elements of a Geometric Algebra are given a geometric interpretation. 
This leads to a focus on operations that are defined on geometrically meaningful 
subsets of the algebra and therefore to the introduction of additional structure 
(and more products between elements) , so that a consistent geometric interpre- 
tation can be maintained on the results of computations. To emphasize this 
difference, we call the elements of the Geometric Algebra muUivectors and we 
call the standard (Clifford) product between elements in the Algebra the geo- 
metric product. The geometric product is denoted by juxtaposition of operands, 
as in AB. 

A summary of the extra features and terminology of Geometric Algebra 
follows. 

1. Blades 

If a nonzero multivector can be written as the geometric product of r 
mutually anticommuting vectors, then it is called an r-blade. The word 
blade refers to an r-blade with the value of r unspecified. Real numbers 
are considered 0-blades and are often called scalars. Vectors are considered 
1-blades. The square of a blade is a scalar. For a vector the square is the 
squared length determined by the metric. Zero is considered an r-blade 
for any value of r, but any nonzero blade is an r-bladc for only one value 
ofr§. 

2. Steps (or Grades) 

A linear combination of r-blades will be called an r-vector, and will be said 
to have step (or grade) r. The space of r-vectors is a linear subspace of 
the entire Clifford Algebra. An arbitrary multivector. A, can be uniquely 
written as J2{^)r where {A)r is the r-vector part of A if r is a nonnegative 
integer and {A)r is zero if r is not a nonnegative integer. 

3. Outer Product 

The outer product of an r-vector A and an s- vector B is defined to be 
{AB)s+r- It is denoted hy AA B and is extended by linearity to arbitrary 
multivectors. The outer product is associative between all multivectors 
and anti-symmetric between vectors. The outer product of two blades is 
a blade (see the appendix). 

4. (Contraction) Inner Product 

The (contraction) inner product of an r-vector A and an s-vector B 
is defined to be {AB)s^r- It is denoted by A\B and is extended by 
linearity to arbitrary multivectors. This inner product differs slightly 
from the inner product of Hestenes|^. It has the useful properties that 
(A A B)\C = A\{B\C) and elA ^ a\ A + a. A A for any multivectors A, 
B, C and any vector a. These details are obvious from considering this 
definition and looking at the proofs in HestenesQ. The inner product is 
explicitly expanded for the vectors a, Ci, C2, ... , and Cfe as 

aj(ci A 02 A ... A Ck) = ^(-l)'^+^(ajcj)ci A ... Ac,. A ... A (1) 
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where the inverted circumflex indicates that the rth vector was omitted 
from the product. The proof in reference carries over to the contraction 
inner product with no modification. The inner product of two blades is a 
blade (see the appendix). 

5. Subspaces 

One of the geometric interpretations common in geometric algebra is to 
use blades to represent subspaces. This works because blades are closely 
related to subspaces. Given a nonzero blade, A, the set {x e 7?." : xA A = 
0} is a subspace of TZ". Similarly given an oriented basis {xi,X2, ...,Xfc} 
for a /c-dimensional subspace A, there is a nonzero blade that corresponds 
to that oriented basis. If fc = then that blade is 1. If fc ^ then that 
blade is Xi A X2 A . . . A Xfc . The identity x G A x A A = is the 

means by which Geometric Algebra subsumes the operation G from set 
theory. Since the square of a blade is a scalar, the inverse of a blade (if it 
has one), is equal to the blade divided by its square, so the inverse is just 
a scalar multiple of the original blade, and hence a blade and its inverse 
represent the same (unoriented) subspace. 

6. Pseudoscalars 

Given an algebra, TZp^q^r, and a set of {p-\-q+r) vectors {vi, Vp+g+^j, the 
outer product vi A ...AVp+q+r is called a pseudoscalar. Often a particular 
nonzero pseudoscalar is singled out. This preferred pseudoscalar serves to 
determine the reference orientation for the vector space and sometimes is 
used to perform duality operations. Usually the preferred pseudoscalar is 
called the pseudoscalar and denoted I. This can be seen as merely a special 
case of the previous section on subspaces, by noticing that every nonzero 
fc-blade is a pseudoscalar for the k dimensional subspace it represents. 

7. Projection Operators 

For a Euclidean metric, blades are closely related to projection opera- 
tors. Given a nonzero blade. A, representing the subspace A, the vector 
^a(x) = (xj A)A^^ is the orthogonal projection of the vector x onto the 
subspace A. As taken fromj|], the identities 



hold for any nonzero blades A and B. Equation (y) implies that a projec- 
tion operator can be decomposed analogously to the way its corresponding 
blades can be factored. 

8. Outermorphism 

Given a linear function, / : T^P^^.f ^ TZ^-i-'^^ there is an extension of 
the function to arbitrary multivectors called the outermorphism of /. It 
is denoted by / and is defined to be the identity when restricted to the 
scalars. The condition /(AA (B + C)) = /(A)A/(B) + /(A)A/(C) is then 
sufficient to define the outermorphism on arbitrary multivectors. This 
extension is well-established [El. 



AB = A A B ^ Pab = Pii + PA 



(2) 



and 



AB = AJB => Pab = Pb - Pa 



(3) 
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A notation convention is adopted to aid the reader in easily making mean- 
ingful distinctions between different multivectors. Lowercase Greek letters are 
reserved for scalars. Lowercase Latin letters are reserved for integers or func- 
tions when not in bold face, while lowercase Latin letters are reserved for vectors 
when in bold face. Bold face is reserved for blades. Lastly, uppercase Latin let- 
ters are used when it is impossible or unnecessary to be more specific about 
the nature of a multivcctor. This notation convention simplifies the reading of 
equations and emphasizes that different geometric interpretations are applied 
to different elements of the algebra. 

3 The Euclidean Metrics 

Since every blade represents an oriented subspace and blades are easy to com- 
pute with, they are a natural candidate for subspace computations. The extra 
scalar degree of freedom allows the future potential for more precise calculations 
with subspaces that attach meaning to the magnitude of a blade. Therefore 
in this paper four blade operations are introduced to correspond to the four 
subspace operations. The blade operations are made first for Euclidean space 
because the relationship between blades and projection operators is strongest 
in Euclidean space. Therefore first a correspondence is made from projection 
operators to blades, then the four blade operations are defined and shown to 
faithfully mirror the subspace operations. 

3.1 The Blade Correspondence 

Here we give an algorithm to construct a blade from its corresponding projection 
operator. The algorithm has an arbitrary scale and orientation inherited from 
an arbitrarily chosen basis, which is the best that can be expected. Let P4 be an 
idempotent linear operator on 7^" such that the image of 7?," is a i-dimensional 
subspace, A. The algorithm constructs a blade that characterizes A as follows. 
First we construct a set of candidate blades, and then show that all the nonzero 
candidate blades represent the subspace A. Finally we show that at least one of 
the candidate blades is, in fact, nonzero. Let Vi, V2, ... , be k f-blades such 
that {Vi, V2, ... ,Vfc} is a basis for the space of vectors. Let Tj = PA(Vi). 
The set {Ti, T2, T^} is our set of candidate blades. By the properties of 
the outermorphism, each Tj is a i-blade. Each is clearly the outer product 
of t vectors and each of these t vectors is in A. If the t vectors are linearly 
dependent then = 0, if not then they form a basis for A so Tj 7^ and Tj 
is exactly the kind of blade to characterize the subspace A. All that remains is 
to show that at least one of the candidate blades is nonzero. Since a blade that 
characterizes A exists we know that it is a linear combination of {Vi , V2, ... 
, Vfe}, so by the linearity of the outermorphism there must be a Vj such that 

Given this correspondence we can translate operations between projection 
operators into operations between blades, except for a loss of the scale and 
orientation. 
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3.2 The Inner Division Operation 



Consider two nonzero blades, A and B, that characterize the subspaces A and 
B respectively. When S is a subspace of A we use the expression A\B to 
denote the quantity B~^A and the call the operation inner division. The inner 
division operation is motivated by the difFerence operation for subspaces. The 
justification requires showing two points, first that under such conditions, B~^A 
is a blade, and second that x A (B~^A) = <^=> x e A\B. 

Let {bi,...,bs} be an orthogonal basis for B. Let {ai,...,ar} be an or- 
thogonal basis for A\B. Clearly {bi, b^, ai, a^.} is an orthogonal basis 
for A. Since {bi,...,b5} is an orthogonal basis for B and since B~^ char- 
acterizes the subspace B it follows that B"-'^ is a nonzero scalar multiple of 
bib2...bs. Similarly A is a nonzero scalar multiple of bib2...bsaia2...ar-. It 
then follows that there exists two nonzero scalars a and (3 such that B~^A = 
Q(bib2...b,5)(bib2...bsaia2...a,.) = /3aia2...cir. Thus B~^A is a blade and it 
characterizes the subspace A\B. 

A quick look at the step of the output reveals that when S is a subspace of 
A then B~^A = B~^JA. Therefore, while it appears that the inner division 
operation is based on the geometric product it is also just as easily based on 
the inner product. It is useful whenever a product between two blades can be 
written as either of two products, because cither definition can be used from 
line to line of a computation, depending on which product gives simplifications 
at that particular moment. An example is the identity, 

A = B(A\B) = B A (A\B) (4) 

which is proved by decomposing the inner division first into the geometric prod- 
uct and then into the inner product and looking at the step of the outcome. 

3.3 The Delta Product 

Consider two nonzero blades, A and B, that characterize the subspaces A and B 
respectively. Let C = AnB and let C be any blade characterizing that subspace. 
When C = {0}, A A B 7^ 0. When C ^ {0}, A A B = 0. In the latter case we 
can define Aj^ = AC and Bj^ = C~-'^B. Since C and C"-'^ both represent C, 
which is a subspace of both A and B, the previous section makes it clear that 
A_L and Bj_ are blades. The intersection of the subspaces characterized by A_l 
and B^ contains only the element zero, so Aj_ A Bj^ 7^ 0. 

By construction AB = Aj_B_l, so the highest step portion of AB is Aj_ABj_, 
and therefore a blade. This motivates a new product for blades which we call the 
delta product. The delta product of two blades, A and B is denoted AAB and 
defined to be the highest step portion of AB. The delta product is motivated 
by the symmetric difference operation for subspaces. The justification requires 
showing that x A (AAB) = x G AAB. 

Let {ai, a^} be an orthogonal basis for A\C. Let {bi, b^} be an orthog- 
onal basis for B\C. Similar arguments as in the previous section demonstrate 
that there exists a nonzero scalar a such that AB = aa.ia.2...SLrhih2...hs, so 
AAB = aaiAa2A...AarAbiAb2A...Ab^. Clearly x A (AAB) = ^ x G AAB. 
Therefore assume x G AAB and we will show that x A (AAB) = 0. By defi- 
nition Ela G ^ and 3h E B such that x = a + b and Vc G C x • c = 0. Let 
ao = a -|- Fc (b) and bg = h + Pc (a) . Clearly ag -|- bo = x + Pc (x) = x. Clearly 
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SiQ E A and agJC = aJC — bJC = 0, so Vc G C, cjao = therefore ao G A\C. 
This means ao A ai A a2 A ... A a^ = so ao A (AAB) = 0. Similarly for bo. 
Therefore x A (AAB) = 0. 

The delta product is very different than the inner or outer product in its 
algebraic properties. The major difference is that AA(B + C) ^ AAB + AAC, 
so the delta product cannot be extended by linearity to arbitrary multivectors. 
The delta product can only be used on blades. More care must be taken with 
implementations of the delta product because even a small change in either A 
or B can cause a change in the step of AAB. 



3.4 The Meet Operation 

Consider two nonzero blades, A and B, that characterize the subspaces A and 
B respectively. The blade correspondence can be used to define a new product 
for blades called the meet and denoted A n B. A n B is defined modulo a scale 
and an orientation as the blade corresponding to the projection operator Pahb, 
where Pahb is defined as: 

„ i^Blx) -PaAb(x) +P(AAB)B-i(x) 

i^AnB(x) = ^ '- (5) 

The justification requires that Pahb be the orthogonal projection onto AClB. 
Let C be a blade characterizing the subspace ACi B, then define A_l = AC 
and B_L = C~^B as above. Now the result follows from a simple application 
of equation (^. First note that Bj^ = B\C, so equation implies that 
B = CB^ = CABx, therefore we have the following identity: 



Pb = Pb^ + Pc (6) 

Similarly A_l A B^ ^ B_l((A_l A Bj.)\Bi) ^ Bx A ((A^ A Bj.)\Bx), 
therefore we have the following identity: 



^^AAB = -Pb-I(aab) + (7) 

Finally (AAB)C^^ = (AAB) A C^^, because every vector in C^^ is or- 
thogonal to every vector in AAB. Since (AAB)B2^ represents a subspace of 
AAB it is just as true that ((AAB)B2i)C-i = ((AAB)B~^) A therefore 
we have the following identity: 

-P(AAB)B-i = ■P(AAB)B];i + -^C-i (8) 

The linear combination of the three projection operators has now been re- 
duced to the linear combination of four projection operators. A quick appeal to 
equation (^ implies that for two blades, if their geometric product is a scalar 
then their projection operators are equal. Now CC~^ is a scalar, and so is 
BJ^(AAB)(AAB)B]|^. Therefore the terms ^'(aab)B"^ ^b"^(aab) 
eel and the terms Pc-i and Pc combine. The result, (equation (||)), then follows 
from equations (^, (^, and (^. 

A small commentary is in order. The first comment is that the calculation 
of the blade correspondence will be simplified by the fact that if the steps of 
A, B, and AAB are r, s, and q respectively, then the step of the meet is 
The second comment is that the blade correspondence is not precise 
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about the scale and orientation of the blade because the blade correspondence 
inherits an arbitrary scale and orientation from an arbitrary basis of blades. 
Since the meet for blades is defined by the blade correspondence, this lack of 
precision is then passed on to the meet for blades, except for the disjoint case. 
The disjoint case occurs when AO B = {0}, and in this case one can choose 
a basis a priori. This is possible because in this case the meet for blades is 
a scalar. Therefore one can choose the scalar '1' for the basis, and then since 
the blade correspondence uses the outermorphism of the projection operator 
and an outermorphism is the identity when restricted to the scalars, the blade 
correspondence gives a determinate answer for the meet, namely it gives '1' back 
again. The third comment is that the meet for blades given here is different from 
previous literature js), which only relates nontrivially to our definition when 
AU B — 72." . When AUB ^ TZ^ the previous literature gives the zero blade as 
the result, while we can also treat that case. The price we pay is linearity. Like 
the delta product, the meet is not linear. 

3.5 The Join Operation 

Previously we noted that (AAB)C = (AAB)AC, in fact (AAB)AC character- 
izes AyjB, therefore using equation (||), we find that Paub = ^'aab+^ahb- We 
have two alternatives. The first alternative is to define the projection operator 
for the join directly as 

n , . ^'(AAB)B-i(x) +PaAb(x) +Pb(x) 

Paub(x) = ^ '- (9) 

and if the steps of A, B, and AAB are r, s, and q respectively, then the 
step of the join is The second alternative is to define the join for blades 

directly in terms of the meet for blades and the inner division operation through 
the equation A U B = A A (B\(A H B)). Just as the meet had a definite 
scale and orientation only in the disjoint case, this definition allows the join 
to inherit the definite scale and orientation from the meet, since in that case 
A U B = A A B. It bears mentioning that this join for blades only agrees with 
the previous literature|^, |j in the disjoint case, but again definitions in the 
previous literature are merely zero when AC\B ^ {0}, so this definition is an 
extension. As with the meet the price we pay is linearity. 

4 The Non-Euclidean Metrics 

A new concept is introduced to allow the subspace operations to be performed 
in any metric. Each of the four operations is then investigated in turn. 

4.1 A LIFT Between Clifford Algebras 

Given two algebras, Ttp^q.r and TZa.b.c, such that a-\-h+c=n, and a 

linear invertible function, /, from the vectors of Ti-p^q.r to TZa,b,c then / is a 
linear invertible map between the two algebras. Call the extended function / a 
LIFT ('linear invertible function' transformation) of 'R'p,q,r to TZa,b,c- 

A LIFT can be used to transfer a problem with subspaces to another algebra, 
preserving incidence relations but allowing the metric to change. Often a LIFT is 
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taken to a Euclidean space. After solving the problem in that space, subspace 
results can be pulled back to the original space. Examples that extend the 
previous results on the subspace operations follow. 

A minor variation of the LIFT is for / to go into an n-dimcnsional subspace 
of a Clifford Algebra over a larger vector space, then the outermorphism is a 
linear invertible map onto a subalgebra. This is especially nice when the outer- 
morphism is an isomorphism between the original algebra and the subalgebra. 
Such a LIFT is called an embedding LIFT (or e-LIFT). This is especially com- 
mon for degenerate algebras TZp,q,r, for which an e-LIFT to TZp+r,q+r,o always 
exists. To see the existence of the e-LIFT, let {pi, Pp, qi, q^, ri, r,.} be 
an orthogonal basis for the vectors in 'R-p,q,r such that pf = 1, q^ = —1, and 
= and let {ei, Gp+j., fi, ■■■,iq+r} be an orthogonal basis for the vectors in 
TZp+r,q+r,o such that ef = 1 and f| = —1. Then let / be a linear function such 
that f{pi) = Bi, f{qj) = fj, and /(r^) = Bp+k + fq+k- Then / is the promised 
isomorphism. 

4.2 The Meet Operation 

If one fixes an arbitrary LIFT, /, from TZp_q_r to TZp+q+r,o,o then the meet, AflB, 
between two blades A and B can be defined by: 

AnB = /-i(/(A)n/(B)) (10) 

The scale and orientation of /(A) fl /(B) are indeterminate except when 
/(A) and /(B) are disjoint, which only happens when A and B are disjoint. 
The LIFT is an outermorphism, so it is the identity on the scalars, so the meet 
has a definite scale and orientation in the disjoint case, and they are the same 
scale and orientation as in the Euclidean space. More importantly, in the disjoint 
case, the scale and orientation are independent of which LIFT, /, was chosen. 
The preservation of the outer product and the scalars makes it clear that this 
meet corresponds to ^ fl S. This means that this definition has a well-defined 
scale and orientation in exactly the cases where the Euclidean definition did, 
and it has an arbitrary scale and orientation in exactly the cases where the 
Euclidean definition did. 

4.3 The Join Operation 

If one fixes fixes an arbitrary LIFT, /, from 7^p,q,r to TZp+q+r,o,o then the join, 
A U B, between two blades A and B can be defined by: 

AUB = /-1(/(A)U/(B)) (11) 

The scale and orientation of /(A) U /(B) are indeterminate except when 
/(A) and /(B) are disjoint, which only happens when A and B are disjoint, 
in which case the join should reduce to the outer product. The LIFT is an 
outermorphism, so it preserves the outer product, so clearly (/(A) A/(B)) = 
/"^(/(AAB)) = AAB. Therefore the join has a definite scale and orientation in 
the disjoint case, and the scale and orientation are independent of which LIFT, 
/, was chosen. The preservation of the outer product and the scalars makes it 
clear that this join corresponds to ^ U B. This means that this definition has 
a well-defined scale and orientation in exactly the cases where the Euclidean 
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definition did, and it has an arbitrary scale and orientation in exactly the cases 
where the Euclidean definition did. 

4.4 The Inner Division Operation 

Consider two nonzero blades, A and B in TZp^q^r, that characterize the sub- 
spaces A and B respectively such that B C A. When exists we can 
calculate A\B = B~^A as usual. Otherwise, we need an e-LIFT, /, from 
T^p,q.r to TZp+r.q+r.o- Let I be the pseudoscalar of TZp+r.q+r.o- Then define 
A\B — /~^(/(A) n (/(B)/)). This meets the definition for the subspace oper- 
ation, but now the scale and orientation has an arbitrary dependence on /. 

4.5 The Delta Product 

Consider two nonzero blades, A and B in TZp^q^r, that characterize the subspaces 
A and B respectively. If AflB has an inverse then the geometric product, AB, is 
nonzero and the highest step portion represents the symmetric difference, AAB 
as usual. Otherwise, we need an e-LIFT, /, from TZp^q^r to TZp+r,q+r,o- Let I 
be the pseudoscalar of TZp+r,q+r,Q- Then define AAB = /"^((/(A) U /(B)) n 
((/(A) n /(B))/)). This meets the definition for the subspace operation, but 
now the scale and orientation has an arbitrary dependence on /. Note that 
the symmetric difference was used to compute the meet and join for Euclidean 
signatures, but the meet and join are used to define the symmetric difference 
for Non-Euclidean signatures. 

5 Linearity and the Meaning of Zero 

The meet and join for blades presented in this paper are not linear, (e.g. in 
general (A-fB)nC7^AnC-KBnC, even when A -f B is a nonzero blade). 
The previous literature[Q ^ have linear results for the meet and join. In our 
notation the meet and join of the previous literature are {AT)\B and AAB 
respectively. 

The linear versions can operate on any multivector (by linear extension), 
but the geometric interpretation of the computation becomes confusing. Also, 
even when the linear versions operate on nonzero blades, they can disagree with 
the subspace operations by giving a result of zero. It is not surprising that the 
subspace operations disagree with the linear operations, because it was exactly 
the preponderance of the answer for many meaningful computations that 
motivated the creation of the new subspace operations of this paper. However 
the opposing versions can be reconciled by pursuing a geometric interpretation 
of the zero blade. 

An interpretation of the zero blade that is consistent with the general en- 
terprise of representing oriented subspaces by blades is to have the zero blade 
represent an indeterminate oriented subspace. At first glance, it appears that 
the zero blade represents the whole space (i.e. {x S TZ"' : x A = 0} = 7?."), 
but this interpretation would imply that the zero blade represents a different 
subspace depending on a particular enveloping pseudoscalar (a property that 
destroys natural subalgcbra and enveloping algebra relationships). But since 
this is implicit in stating that the zero blade can represent any subspace this 
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is actually support for the interpretation proposed here. Furthermore, since 
the hope of representing subspaces by blades is to eventually be able to deal 
with uncertainty in geometrical computations, the scale factor would naturally 
be used to represent how well-determined the blade is. This implies that the 
zero blade represents a completely undetermined subspace. Lastly the linear 
versions give the zero blade as the result if either of the input blades is the zero 
blade, this adsorbing property is consistent with the indeterminacy the zero 
blade represents. 

When two oriented subspaces do not span the pseudoscalar, the linear meet 
gives the result of the zero blade because the orientation of the meet cannot 
be determined. The linear meet is a fully functional quantitative operation, 
which can give a quantitative meet, but only if given a quantitative join first 
(in the role of the pseudoscalar) . Similarly when two oriented subspaces have a 
nontrivial intersection, the linear join gives the result of the zero blade because 
the join has an undetermined orientation. 

With this interpretation for the zero blade, the linear meet and join can 
be compared to the versions presented in this paper. The linearity can be 
an advantage for implementation for some applications, and if that advantage 
outweighs the costs of getting the zero blade as a result, then an educated 
decision to implement the linear versions can be made for that application. 
The interpretation for the zero blade can also be used to extend the subspace 
operations for nonzero blades to the zero blade by declaring the inner division, 
the delta product, the meet, and the join to be zero if cither input of the two 
input blades is the zero blade. 

6 Conclusion 

Blades can represent subspaces, and in applications we need to perform oper- 
ations on subspaces. Therefore we would naturally want operations on blades 
that mirror the results of computations that we would have liked to perform on 
subspaces. The four subspace operations (inner division, delta product, meet, 
and join) are the first four operations from nonzero blades to nonzero blades. 
The hope is that these operations can contribute to quantitative computations 
with oriented subspaces. 

The four operations are on equal footing because the delta product and the 
inner division are used to define the meet and join in Euclidean signatures, while 
the meet and join are used to define the delta product and inner division in Non- 
Euclidean signatures. This indicates that they are all fundamental (though not 
as fundamental as the geometric product) and that they are all interconnected. 

Standard concepts in geometric algebra needed to be augmented because the 
meet and join for blades cannot be defined|^ with an orientation due to fun- 
damental geometric problems. This fundamental problem was solved by using 
projection operators to represent unoriented subspaces. Within this solution 
the delta product helps to compute the meet and join for blades. 

The price to be paid for this augmentation is that the new blade operations 
are not linear and cannot be extended to arbitrary multivectors. The authors 
believe that their non-linearity might be the reason that the operations have 
not been used previously. It is only by sacrificing linearity, and thus losing 
applicability to arbitrary multivectors, that one can solve the meet and join for 
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blades. 

The four subspace operations are tools intended for general use in applica- 
tions, however no applied examples are included in this paper. Readers looking 
for examples of the inner division and the delta product need look no farther 
than the proof of the meet for Euclidean metrics, and readers looking for exam- 
ples of the meet and join need look no farther than the inner division and delta 
product for Non-Euclidean metrics. Beyond the four subspace operations, this 
paper utilizes another tool of general applicability. This tool is the LIFT ('linear 
invertible function' transformation). It is an invertible map between algebras of 
different signatures. This tool can be used to advantageously transform prob- 
lems that are independent of signature to whichever signature is most helpful 
at any particular moment. This tool is demonstrated in the paper by extending 
the results of the meet and join from the Euclidean case to the Non-Euclidean 
case (even to degenerate signatures). 

Lastly some loose ends are resolved. A geometric interpretation is given to 
the zero blade that explains the results of the previous literature and facilitates 
the educated choice between different versions of the subspace operations. The 
final loose end is resolved by the appendix, which includes proofs to demonstrate 
that the inner and outer products go from blades to blades (even in degenerate 
signatures). 

7 Appendix 

In this Appendix we prove that the outer product of two blades is a blade and 
that the inner product of two blades is a blade. To our surprise, this does 
not appear to have been proved before, but is less trivial than may have been 
assumed when degenerate algebras are considered. 

7.1 The Outer Product 

By associativity of the outer product it suffices to show that the outer product of 
a vector and an r-blade is a blade. The result is trivial because the (r-l-l)-vector 
determines an (r-l-l)-dimensional subspace and that subspace has an orthogonal 
basis. However a more constructive proof is desirable, to assist in the proof for 
the inner product and to see how such a factorization can be made. The proof 
proceeds by induction on r, the base case is r = 1, or the outer product of two 
vectors is a blade. 

Let Vi and V2 be two vectors. Since i(vi — V2)(vi -I-V2) = Vi Av2 -I- , 
this gives a factorization of Vi A V2 when •v\ = . If ^ v| then either or 
V2 is not equal to 0. Since vi A V2 = — V2 A vi, we may assume without loss of 
generality that ^ 0. Then we note that vj~^J(vi A V2) is a vector and that 
Vi(vj"^J (vi A V2)) = vi A V2, so this gives a factorization of Vi A V2. Therefore 
in both cases the outer product of two vectors can be factored. 

Assume that the outer product of r vectors is a blade. Let {a, Ci, 0^} be 
r+1 vectors. The inductive step has three cases. 

1. The first case is when ajc^ = for each i . Then by equation (|^), a A Ci A 
...Ac,. = a(ci A... ACr). By the inductive hypothesis, Ci A... Ac^ is a blade, 
so there exist r anticommuting vectors, {ai, ...,ar}, such that aia2...ar = 
Ci A ... A Cj.. Each a^ is a linear combination of the set {ci, Cr}, so each 
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Hi anticommutes with a, so aai...ar is a factorization of a A Ci A ... A c^., 
therefore a A Ci A ... A is an (r+l)-blade. 

2. The next case is when a^ 7^ 0. Let = a~^(aACi). Then aAci A ... Ac^ = 
a A bi A ... A hr- Now we have guaranteed that ajb^ =0 for each i, so the 
previous case resolves the factorization. 

3. The last case is when there exists a k such that ajc^ 7^ and a^ = 0. 
Without loss of generality we may assume k = 1 since the order of the 
vectors only affects the sign of the outcome. By the base case, a A Ci is 
a 2-blade. Direct computation shows that the square of a A Ci is (ajci)^, 
hence nonzero, therefore there exist two invertible vectors c and di such 
that a A ci = cdi. Let di = (cdi)~^((cdi) A Ci) for i G {2,...,r}. Then 
a A ci A ... A c,. = c Adi Ad2 A ... Adr. Now since 7^ the previous case 
resolves the factorization. 

7.2 The (Contraction) Inner Product 

Let A = ai A a2 A ... A a.,„ and B be two blades in TZp.q.r- By the properties 
of the inner product, AJB — aij (a2j (...J (a^JB)...)), so it suffices to show that 
aj C is a blade for any vector a and any blade C. If C is a scalar then aj C = 
so it is a blade. Otherwise let C = Ci A C2 A ... A c^. The nonzero blade C 
characterizes a subspace with signature {s,t,u) and hence a subalgebra TZs,t,u- 
Using equation (|l|), it is clear that aJC resides in TZs^t,u- Let / be a LIFT from 
T^s,t,u to TZkfifi- Since f{X) is a blade if and only if X is a blade then to show 
that aJC is a blade it suffices to show that /(aJC) is a blade. From equation 
(P, it follows that /(aJC) is the sum of (fc-l)-blades. Let I be a nonzero 
pseudoscalar for TZk,n.o- Let {Hi,...,H„i} be m (fc-l)-blades in TZk.oft- Then 
Hi + ... + H,„ = (Hi + ... + H™)I-il = vl, where v = Hil-^ + ... + 'h„J-i = 
HiJI^-'^ + ... + HmJI^-'^ is a vector. If v = then the sum, Hi + ... + Hm = vl, is 
zero and trivially a blade. If not then v 7^ and therefore v has an inverse and 
by the details of the proof for the outer product, it is clear that I can be factored 
under the geometric product with v^^ as a factor. Then Hi + ... + H,„ = vl is 
clearly a blade. 
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